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ABSTRACT 

We present 12 CO(J=1^0) observations of the high-redshift quasi-stellar objects (QSOs) BR 1202- 
0725 (z = 4.69), PSS J2322+1944 [z = 4.12), and APM 08279+5255 (z = 3.91) using the NRAO 
Green Bank Telescope (GBT) and the MPIfR Effelsberg 100 m telescope. We detect, for the first 
time, the CO ground-level transition in BR 1202-0725. For PSS J2322+1944 and APM 08279+5255, 
our observations result in line fluxes that are consistent with previous NRAO Very Large Array 
(VLA) observations, but they reveal the full line profiles. We report a typical lensing-corrected 
velocity-integrated intrinsic 12 CO(J=1^0) line luminosity of L' co = 5 x 10 10 K krus _1 pc 2 and a 
typical total H2 mass of = 4 x 10 10 Mq for the sources in our sample. The CO/FIR luminosity 

ratios of these high-z sources follow the same trend as seen for low-z galaxies, leading to a combined 
solution of log(LFm) = (1-39 ± 0.05) x log(Lco) — 1-76. It has previously been suggested that the 
molecular gas reservoirs in some quasar host galaxies may exhibit luminous, extended 12 CO( J=l— >0) 
components that are not observed in the higher - J CO transitions. Utilizing the line profiles and the 
total intensities of our observations and large velocity gradient (LVG) models based on previous results 
for higher— J CO transitions, we derive that emission from all CO transitions is described well by a 
single gas component where all molecular gas is concentrated in a compact nuclear region. Thus, our 
observations and models show no indication of a luminous extended, low surface brightness molecular 
gas component in any of the high-redshift QSOs in our sample. If such extended components exist, 
their contribution to the overall luminosity is limited to at most 30%. 

Subject headings: galaxies: ISM — cosmology: observations — galaxies: active — galaxies: starburst 
— galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

Understanding when and how galaxies form is one of 
the primary objectives in both observational and theo- 
retical astrophysics. The mere fact that active galaxies, 
such as radio galaxies, QSOs, and emission line galaxies 
are observed up to redshifts of z = 6.6 (e.g., Fan et al. 
2001; Rhoads et al. 2001; Hu et al. 2002; Kodaira et al. 
2003; Taniguchi et al. 2005), less than one Giga-year af- 
ter recombination, implies that galactic scale (~10kpc), 
gravitationally bound structures exist at this early epoch. 
The study of the masses and dynamical state of these 
young systems serves as a direct constraint to the mod- 
els describing the growth of large scale structures since 
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the epoch of recombination. 

Studies of the molecular and dusty interstellar medium 
(ISM) in these galaxies are of fundamental importance, 
since it is this medium out of which stars form; accurate 
determination of the molecular gas mass may therefore 
serve as an indicator of the evolutionary state of a galaxy. 
The detection of carbon monoxide (CO) is also a strong 
confirmation that star formation is going on in some of 
the highest redshift systems. In fact, the combination of 
molecular gas and dust detections with large far-infrared 
(FIR) luminosities provides the strongest evidence that a 
significant fraction of high-z galaxies is undergoing star- 
bursts at prodigious rates (> 10 3 M Q yr _1 ), consistent 
with the formation of a large elliptical galaxy on a dy- 
namical timescale of ~10 7 -10 8 yr. 

Over the past decade, more than 30 galaxies at z > 2 
have been detected in CO emission (Solomon & Vanden 
Bout 2005) out to a redshift of z = 6.42 (Walter et al. 
2003, 2004; Bertoldi et al. 2003), confirming the pres- 
ence of intense starbursts in numerous high-z galaxies. 
As most of these observations were obtained using mil- 
limeter interferometers, these detections were typically 
achieved by observing high-J CO(J — ► J — 1) transi- 
tions (J > 3). Although these high-J lines exhibit in 
general higher peak flux densities than the ground-state 
(J = 1) transition, it is possible that the higher order 
transitions are biased to the excited gas close to a central 
starburst and do not necessarily trace the entire molec- 
ular gas reservoir seen in 12 CO(J=1^0). Also, the con- 
version factor ('a') to derive molecular (H2) masses from 
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measured CO luminosities has mostly been estimated for 
the 12 CO(J=1^0) line (e.g., Downes & Solomon 1998; 
Weiss et al. 2001). Observing 12 CO(J=1^0) has the 
additional advantage that properties of the highly red- 
shifted sources can be directly compared to the molec- 
ular gas properties of nearby (starburst) galaxies which 
are predominantly mapped in the 12 CO( J=1^0) transi- 
tion. 

These are main motivations for observing the 
12 CO( J=l— >0) ground-state transition. However, due to 
the faintness of the line and the bandwidth limitations of 
current radio telescopes, such high- z 12 CO( J=l— >0) ob- 
servations only exist for two QSOs and two radio galax- 
ies to date. All of these observations have been ob- 
tained with radio interferometers operating at centime- 
ter wavelengths - the NRAO Very Large Array (VLA; 
Papadopoulos et al. 2001; Carilli et al. 2002a; Greve et 
al. 2004) and the Australia Telescope Compact Array 
(ATCA; Klamcr et al. 2005). In particular, due to the 
bandwidth limitations of the VLA, obtaining better con- 
straints on the spectral line shape and total flux of the 
12 CO(J=1^0) transition is desirable even for already 
detected sources. Today's largest single dish telescopes, 
such as the NRAO Green Bank Telescope (GBT) 10 and 
the MPIfR Effelsbcrg telescope 11 , can eliminate some of 
those issues due to their larger spectral band widths. 

Here we report on first observations of 12 CO( J=l— >0) 
in three CO-bright high-z QSOs, which also exhibit 
ultra-luminous IR emission, using the GBT and the Ef- 
felsbcrg telescope. In Section 2, we describe our observa- 
tions. Section 3 summarizes our results on the individual 
objects (BR 1202-0725 at z = 4.69, PSS J2322+1944 at 
z = 4.12, and APM 08279+5255 at z = 3.91). Section 4 
provides an analysis and discussion, and Section 5 closes 
with a summary of our results. We assume a standard 
A cosmology throughout, with H = 71 kms -1 Mpc -1 , 
n M = 0.27,' and Q A = 0.73 (Spergel et al. 2003). 

2. OBSERVATIONS 
2.1. GBT 

Observations of the three targets were carried out with 
the GBT during 12 observing runs between October 2004 
and April 2005 with a total observing time of 89.5 h. 
The total/on-source observing times were 31.5/20 h for 
BR 1202-0725, 23/15h for PSS J2322+1944, and 35/22h 
for APM 08279+5255. During all runs, 3C147 and 3C286 
were used as primary/flux calibrators. For spectral line 
calibration, we observed IRC+10216, DR21, OrionIRc2, 
Orion HC, and W30HH 2 0. As secondary/pointing 
calibrators, 3C273, J1256-0547 (for BR 1202-0725), 
J2253+1608 (for PSS J2322+1944), J0753+538, and 
J0824+5552 (for APM 08279+5255) were utilized. The 
pointing accuracy, determined by continuum cross scans 
of nearby sources, was typically ~3" (reaching <1" 
under the best conditions). We estimate the calibra- 
tion to be accurate to 10-15%. As the 12 CO(J=1^0) 
transition at 115.2712 GHz is rcdshiftcd into the K 

10 The Green Bank Telescope is a facility of the National Radio 
Astronomy Observatory, operated by Associated Universities, Inc., 
under a cooperative agreement with the National Science Founda- 
tion. 

11 The Effelsberg telescope is a facility of the Max-Planck- 
Gesellschaft (MPG), operated by the Max-Planck-Institut fur Ra- 
dioastronomie (MPIfR). 



band for all three targets (BR 1202-0725: ^ bs = 
20.2450 GHz, PSS J2322+1944: ^ obs = 22.5258 GHz, 
APM 08279+5255: v ohs = 23.4663 GHz), the dual-beam, 
dual-polarization 18-26 GHz receiver was used for all ob- 
servations. The beam size of the GBT at our observing 
frequencies is ~32-36" (~225-250kpc at z ~ 4), i.e., 
much larger than our targets. The two beams have a 
fixed separation of 178.8" in azimuth direction. Two dif- 
ferent spectrometer setups were used; half of the observ- 
ing runs were executed in the first mode, and the other 
half of the runs in the second mode. The first mode 
features two simultaneous intermediate frequencies (IFs) 
with a bandwidth of 800 MHz (-10,000-12,000 km s" 1 
at our observing frequencies) and 2048 channels each, re- 
sulting in a spectral resolution of 391kHz (—5 kins" 1 ). 
The 12 CO( J=l— >0) line was always centered in the first 
800 MHz IF. The second mode has one IF with 200 MHz 
(~2, 500-3, 000 kms -1 ) and 16384 channels, resulting in 
a spectral resolution of 12 kHz (-0.15 kms^ 1 ). BR1202- 
0725 and APM 08279+5255 were observed with both se- 
tups, while PSS J2322+1944 was only observed during 
runs that used the second setup. The spectra taken in 
each mode were examined separately and then combined. 
The ON-OFF position switching mode was used, i.e., the 
target was observed alternatively with the two telescope 
beams, and the off-source beam was always monitoring 
the sky background in parallel. The beam switching fre- 
quency was once every 60-120 seconds depending on the 
observing run. The weather was excellent for the win- 
ter nights with typical zenith system temperatures of 
T sys = 22 - 35 K on a T£ scale! On April 19 and 20, 
2005, the system temperatures were significantly higher 
(typically T sys = 50 - 60 K) . 

For data reduction, the aips++ 12 package and 
the new GBT IDL software 13 were used, provid- 
ing consistent results. The spectra of BR 1202-0725, 
PSSJ2322+1944 and APM 08279+5255 were binned 
to 2.00MHz (SOkms- 1 ), 1.80 MHz (24 kms" 1 ), and 
5.86 MHz (75 kms -1 ), reaching rms noise values of 
~75^iJy, ~140/xJy, and ~65^iJy. A linear baseline was 
subtracted from the PSS J2322+1944 spectrum to re- 
move continuum fluxes and atmospheric/instrumental ef- 
fects. For BR 1202-0725 and APM 08279+5255, polyno- 
mials of order 2 have been used to remove a very wide 
(3>100MHz, i.e., much broader than the width of the CO 
line) 'bending' of the baselines. Only the frequency range 
that was used to define the spectral baselines is used for 
the final GBT spectra. APM 08279+5255 was observed 
using both the 200 MHz and 800 MHz bandwidth setups. 
While the 12 CO(J=1^0) line was clearly detected with 
both setups, the spectrum obtained with the narrow- 
band high-resolution 200 MHz setup shows significant 
baseline problems, as the continuum level is significantly 
different on the red and blue sides of the line. Therefore, 
we only use the results obtained with the wide 800 MHz 
bandwidth setup (corresponding to 15 h on-source, and 
covering a velocity range of ~8,000 kms -1 ) for the fi- 
nal spectrum. In addition to simple ON-OFF combina- 
tions and subtraction of low-order polynomial spectral 
baselines, we also followed the scheme proposed by Van- 
den Bout et al. (2004) for reconstructing the temporal 

12 aips2.nrao.edu 

13 gbtidl.sourceforge.net 
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baseline variations, however found that this did not im- 
prove our final results. The reduced spectra were read 
into GILDAS/CLASS 14 to write out the final combined 
spectrum tables. The final spectra are shown in the top 
panels of Figs. H and [H and in Fig. El We note that 
the spectrum of BR 1202-0725 is limited to the central 
— 150 MHz due to the fact that about half of the data 
were taken with the narrower 200 MHz setting. 

2.2. Effelsberg 

Observations were carried out towards BR 1202-0725 
and PSS J2322+1944 in January and February 2003 with 
a total observing time of — 40h (20 h per source, cor- 
responding to 8h on-source each). At 40" beam size, 
the pointing accuracy, as determined by continuum cross 
scans of nearby sources, was better than 10", with typi- 
cal values of ~5". Calibrations for the gain as well as the 
variation of the atmospheric opacity and zenith distance 
were obtained from observations of 3C286 and NGC 7027 
(see Ott et al. 1994 for reference fluxes), leading to a to- 
tal formal calibration uncertainty of ~15%. We observed 
in beam switching mode utilizing a rotating horn with a 
beam throw of 2' and a switching frequency of — 1 Hz. 
We have used a dual polarization HEMT receiver for 
all observations. The autocorrelator backend was split 
into eight bands of 160 MHz bandwidth and 128 chan- 
nels each that could individually be shifted in frequency 
by up to ±250 MHz relative to the recessional velocity of 
the targets. The final spectra cover a velocity range of 
~6,000 kms -1 with channel spacings of ~16kms _1 . We 
achieved single channel system temperatures of 65 K and 
85 K on a T£ scale. After combination of both orthogo- 
nal linear polarizations, this leads to a T sys of 47 K and 
60 K. 

The GILDAS / CLASS package was used for data reduc- 
tion. All spectral baselines are of good quality and only 
first order polynomial baselines had to be subtracted. 
The final spectra of BR 1202-0725 and PSS J2322+1944, 
are shown in the middle panels of Figs. ^ and [3 
These spectra were binned to 5.97MHz (88kms _1 ) 
and 3.60 MHz (48kms -1 ), reaching rms noise values of 
— 150/iJy and ~380/iJy. We note that the latter do 
not show the full ranges used for spectral baseline fit- 
ting. Those ranges (-380 MHz for BR 1202-0725, and 
-180 MHz for PSS J2322+1944) are illustrated by the 
insets in the same figures. 

2.3. Combined spectra 

As BR 1202-0725 and PSS J2322+1944 were observed 
with both telescopes, we also created combined spec- 
tra of both results to increase the signal-to-noise ratio 
and average out part of the calibrational uncertainties. 
The GILDAS / CLASS package was used to reprocess and 
combine the final spectra of both telescopes. The final 
unbinned GBT and Effelsberg spectra were regridded to 
a common velocity resolution and baseline-subtracted 
before combination. In the combination, spectra were 
weighted with their respective rms. The combined spec- 
tra of BR 1202-0725 and PSS J2322+1944, as shown 
in the bottom panels of Figs. ^ and |3 were binned 
to a resolution of 2.00 MHz (SOkms^ 1 ) and 1.80 MHz 

14 www.iram.fr/IRAMFR/GILDAS 



(24kms ), reaching rms noise values of — 70/xJy and 
~125/iJy. 

3. RESULTS 

The sources in our study are the three CO-brightest 
high-redshift QSOs that can currently be observed in the 
12 CO( J=l— >0) transition. We obtained detections for all 
our targets. CO line luminosities L' co (in K kms _1 pc 2 ) 
were derived utilizing: 

L' co = 3.25 x 10 7 x I x u~& x Dlx (1 + z)~ 3 (1) 

where I is the velocity-integrated 12 CO(J=1^0) line 
flux in Jy kms -1 , Z?l is the luminosity distance in Mpc, 
and fobs is the observed frequency in GHz (Solomon et 
al. 1992). For the systems with known lensing magnifica- 
tion factor fi^° , L' co has to be divided by that factor in 
order to get the intrinsic CO luminosity of the discussed 
target. A conversion factor a = 0.8 M Q /K kms _1 pc 2 
to convert ^co(i-o) *° -^gasCHs) is assumed throughout, 
as applicable for local ultra-luminous infrared galaxies 
(ULIRGs)/starbursts (Downes & Solomon 1998). All ob- 
servational results are summarized in Tabled Derived 
CO luminosities and gas masses as well as FIR luminosi- 
ties from the literature are given in Tabled 

3.1. BR 1202-0725 

3.1.1. Previous results 

BR 1202-0725 (z = 4.69) was detected in multi- 
ple CO transitions before (see Table |3J), but not in 
12 CO(J=1^0). This optically bright, radio-quiet QSO 
has the curious property that the optical QSO is a sin- 
gle source, but the mm continuum and CO line obser- 
vations show a double source with a separation of about 
4" (Omont et al. 1996; Guilloteau et al. 1999; Carilli 
et al. 2002b). This double morphology may indicate a 
pair of interacting objects separated by only 28 kpc (Yun 
et al. 2000; Carilli et al. 2002b). An alternate expla- 
nation would be a double starburst system composed of 
a QSO (southern source) and a dust-obscured, Ly-a - 
emitting companion (northern source, Hu et al. 1996), 
which is ionized by the strong QSO. Recently, Klamer et 
al. (2004) suggested that the nature of the double source 
might be due to jet-induced star-formation, where the 
northern component corresponds to a radio hot spot. It 
has been discussed that the quasar activity in this sys- 
tem may be triggered by gravitational interaction with 
the companion. Using typical conversion factors for a 
ULIRG (Downes & Solomon 1998), the total molecular 
gas mass derived from the CO luminosity exceeds the dy- 
namical mass of the system (Carilli et al. 2002b), and as 
a remedy for this inconsistence, gravitational lensing has 
been suggested. However, the same authors argue that 
masses and velocity widths of the components are very 
different, therefore it appears unlikely that the multiple 
components are different lensed images of one source. 

3.1.2. New observations 

The final 12 CO(J=1^0) spectra of BR 1202-0725 are 
shown in Fig. n(^°P panel: GBT, middle panel: Effels- 
berg, bottom panel: combined spectrum). Gaussian fit- 
ting to the line profile of the combined spectrum results 
in a peak flux density of 0.34 ± 0.03 mJy and a FWHM 
of 333 ± 30 kms" 1 . The integrated 12 CO(J=1^0) line 
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flux is 0.120 ± 0.010 Jy kms" 1 . This agrees well with 
the extrapolated value of 0.123 ± 0.013 Jy kms" 1 de- 
rived from the 12 CO(J=2^1) flux (Carilli et al. 2002b) 
under assumption of fully thermalized and optically thick 
CO emission (see also discussion in Sect. 4.1). The width 
of the 12 CO( J=l— >0) line is consistent within the error 
bars with an average of the higher-J transitions in the 
literature (~ 290 kms -1 , see Table [3J ■ The derived red- 
shift of 4.6949 ± 0.0003 is in good agreement with the 
12 CO(J=5^4) redshift (Omont et al. 1996, northern and 
southern components combined). The structure of this 
source is unresolved at the resolution of our measure- 
ments; all line parameters are in good agreement with 
the higher-J lines adding up both components. 

The derived CO line luminosity of 1.0 x 10 11 K 
kms^ 1 pc 2 results in an H 2 gas mass of 8.1 x 10 10 M Q 
(see also Tab.EJl. 

3.2. PSS J2322+19U 

3.2.1. Previous results 

The z = 4.12 QSO PSS J2322+1944 is an IR-luminous 
high-redshift source (see Table that is known to ex- 
hibit strong CO line emission in multiple transitions (see 
Table |3J). Optical imaging and spectroscopy reveal a 
double source structure, and the two components are 
separated by about 1.5". The spectra of both peaks 
are essentially identical, consistent with strong gravita- 
tional lensing (optical magnification factor /z^ 5 * = 3.5) 
by an intervening foreground galaxy (Carilli et al. 2003). 
In high-resolution VLA images, Carilli et al. (2003) 
find a molecular Einstein Ring with a diameter of 1.5" 
in 12 CO(J=2— >1) line emission, which can be mod- 
eled as a circumnuclear star-forming disk with a ra- 
dius of 2.2 kpc (CO magnification factor zt^ ~ 2.5). 
The derived intrinsic star-formation rate (SFR) is of 
order gOOMgyr" 1 . PSS J2322+1944 exhibits strong, 
non-thermal (synchrotron) radio continuum emission at 
1.4 GHz, and the rest-frame radio-to-IR spectral energy 
distribution (SED) resembles that of local nuclear star- 
burst galaxies like M82 (Cox et al. 2002). This QSO is 
the fourth high-z CO emitter to be detected in [C i] emis- 
sion (Pety et al. 2004), providing additional evidence for 
the presence of active star-formation in the host galaxy. 

3.2.2. New observations 

The final 12 CO(J=1^0) spectra of PSS J2322+ 1944 
are shown in Fig. {top panel: GBT, middle panel: 
Effelsberg, bottom panel: combined spectrum). From 
Gaussian fitting to the combined spectrum, the peak 
line flux density is found to be 0.77 ± 0.07 mJy, the 
line FWHM is 190 ± 14 kms" 1 , and the integrated line 
flux is 0.155 ± 0.013Jy kms" 1 . These values are in 
good agreement with the 12 CO(J=1^0) detection of 
Carilli et al. (2002a; S v = 0.89 ± 0.22 mJy, AVfwhm = 
200 ± 70 kms" 1 , / = 0.19 ± 0.08 Jy kms" 1 ), although 
results for the higher-J CO transitions indicate a larger 
linewidth (AI/ FWHM > 250 kms" 1 , see Table EJ). The 
Gaussian fit gives a redshift of 4.1179 ± 0.0002. 

The derived lensing-corrected CO line luminosity of 
4.2 x 10 10 K kms _1 pc 2 results in an H2 gas mass of 3.4 x 
1O 1O M (see also Tab.EJ. 

3.3. APM 08279+ 5255 



3.3.1. Previous results 

APM 08279+5255 is a strongly lensed, radio-quiet 
broad absorption line (BAL) QSO at z — 3.91. Gravita- 
tional lens models of the QSO continuum source suggest 
magnification by a factor of /i 1 ^* ~ 100, and the image 
breaks up into three components with a maximum sep- 
aration of 0.4" (Ledoux et al. 1998; Ibata et al. 1999; 
Egami et al. 2000). APM 08279+5255 has been detected 
in the mm and sub-mm dust continuum, revealing an ap- 
parent bolometric luminosity of ~ 5 x 10 15 L Q (Lewis et 
al. 1998). A multi-transition CO study (see Table EJ) in 
combination with detailed lens models appear to reveal a 
spatially extended structure on a scale of at least 400 pc 
(Lewis et al. 2002), which is gravitationally magnified by 
a factor of zt£° = 7. The strength of the 12 CO(J=9^8) 
emission indicates the presence of hot dense gas with a 
kinetic temperature of approximately 200 K (Downes et 
al. 1999). 

From VLA imaging of the CO ground-state transi- 
tion at a linear resolution of 2.25", Papadopoulos et 
al. (2001) report the detection of an extended, low- 
excitation molecular gas reservoir around the compact 
nucleus which extends over a scale of 7" (^30 kpc). The 
integrated brightness temperature of this extended do- 
main appears to be of the same order of magnitude as 
that of the nuclear region. This extended reservoir would 
be well within our 32" GBT beam. Papadopoulos et al. 
(2001) do not derive the total flux in the extended reser- 
voir. For the central ~1" (corresponding to ^7.2 kpc 
at the source redshift), which they call "the nuclear 
12 CO(J=1^0) emission", they find an integrated flux 
of 0.150 ± 0.045 Jy kms" 1 . The VLA bandpass used by 
Papadopoulos et al. (2001) has an effective bandwidth of 
~45MHz, or ~575kms -1 at the 12 CO(J=1-^0) line fre- 
quency. Therefore, assuming that the peak in their Fig- 
ure lc is 6.5cr = 260 /it Jy beam" 1 (a — 40/iJy beam" 1 ), 
we consistently derive an integrated flux over the peak of 
0.15Jy kms -1 beam" 1 . However, Papadopoulos et al. 
(2001) suggest that this peak of 'nuclear emission' sits 
on a broad, extended plateau (see their Figure lc). If we 
assume that the extended reservoir component is traced 
by their 2er (80/iJy beam -1 ) contour, it has a width 
of ^8"x3", or ^4.7 beam areas at the given resolution 
of 2.25" x 2.25". This corresponds to a flux of 0.22 Jy 
kms" 1 for the extended emission. The same estimate 
for their 3<r (120 /iJy beam" 1 ) contour gives a width of 
~7"x2.25", or ^3.1 beam areas. This again corresponds 
to a flux of 0.22 Jy kms" 1 (note that both these esti- 
mates are lower limits for the extended emission). The 
total flux in the 12 CO( J=l— >0) map of Papadopoulos et 
al. (2001) is therefore estimated to be >0.37Jy kms" 1 . 
The same authors also suggest that this reservoir breaks 
up into multiple components at higher resolution; in 
their 12 CO( J=2— >1) observations at 0.5" resolution, they 
claim to find two emitting regions 2-3" distant from the 
central region - if real, these could be companion galax- 
ies which are not individually resolved in the extended 
12 CO(J=1^0) reservoir. 

3.3.2. New observations 

The final 12 CO(J=1^0) spectrum of 
APM 08279+5255 is shown in Fig. We derive a 
12 CO(J=1^0) peak flux of 0.26 ± 0.04 mJy and an 



CO(l-O) in z>4 QSOs 



5 



FWHM of 556 ± 55 kms" 1 from our Gaussian fit. This 
results in an integrated line flux of 0.152 ± 0.020 Jy 
kms , which is in good agreement with the value of 
0.150 ± 0.045 Jy kms' 1 found by Papadopoulos et al. 
(2001) for the central ~1" ("the nuclear 12 CO(J=1^0) 
emission"). However, our result is inconsistent with the 
integrated flux of 0.37 Jy kms" 1 given by our estimate 
of their full 12 CO(J=1^0) reservoir (see previous 
subsection). Such a high flux is clearly ruled out by 
our GBT observations. We thus find no evidence for 
a luminous extended halo, which would be well within 
our 32" beam, and cannot confirm the existence of 
bright companion galaxies. The derived FWHM velocity 
width of our 12 CO(J=1^0) detection is in agreement 
with single-dish observations of the 12 CO(J=6^5), 
12 CO(J=10^9), and 12 CO( J=ll->10) transitions 
obtained with the IRAM 30m telescope (^500 kms , 
Weiss et al. 2006, in prep.) and IRAM Plateau de 
Bure interferometer observations of the 12 CO(J=4^3) 
transition (480 ± 35 kms -1 , Downcs ct al. 1999, see 
Table©. Our derived 12 CO(J=1^0) redshift of 3.9122 
± 0.0007 is in good agreement with previous results 
(3.9114 ± 0.0003 for the 4^3 transition). The derived 
lensing-corrected CO line luminosity of 1.4 x 10 10 K 
kms _1 pc 2 results in an H2 gas mass of 1.1 x 10 10 M© 
(see also Tab.|2}. 

4. ANALYSIS AND DISCUSSION 

4.1. LVG modeling 

To investigate how much of the 12 CO(J=1^0) emis- 
sion in our target QSOs is associated with the molec- 
ular gas reservoirs detected in the higher-J CO transi- 
tions, we have used spherical, one-component large ve- 
locity gradient models (LVG). All these LVG calcula- 
tions use the collision rates from Flower & Pineau des 
Forets (2001) with an ortho/para H2 ratio of 3 and 
a CO abundance per velocity gradient of [CO]/Vu = 

1 x 10" 5 pc(kms- 1 )- 1 (e.g. Weiss et al. 2005). Models 
were fitted to those lines above the 12 CO(J=1^0) tran- 
sition listed in Table 01 for each source. The turnover of 
the CO line SED (and therefore the slope beyond the 
turnover) is not well determined. Thus, a large degen- 
eracy exists between the kinetic gas temperature 7km 
and density /O gas (H2), the two main free parameters in 
our study. As an example, Figure 0] shows data of all 
transitions and three representative models for BR 1202- 
0725: Model 1 (solid line) assumes Tk; n = 60 K and 
p g as(H2) = 10 41 cm~ 3 , and gives the overall best fit 
to all transitions. Models 2 (dashed line) and 3 (dot- 
ted line) are shown as a representation of the parameter 
space allowed by the data within the error bars. Model 

2 assumes T kin = 120 K and p gas (H 2 ) = 10 37 cm~ 3 , 
while model 3 assumes Tki n = 30 K and p gas (H2) = 
10 4 6 cm- 3 . The LVG predicted 12 CO(J=1->0) flux of 
the different models based on the J > 1 CO transi- 
tions is fairly well constrained by the solutions. Most 
of our calculated models suggest that the CO emission 
is close to thermalized up to the 12 CO(J=4^3) tran- 
sition and optically thick (r > 5), which implies that 
the LVG predicted 12 CO(J=1-^0) line fluxes are sim- 
ilar to those we would derive by assuming a v 2 scal- 
ing of the line flux densities from the mid-J CO tran- 
sitions. The predicted LVG integrated flux ranges of 



the 12 CO(J=1^0) transition are 0.10-0.12Jy kms' 1 
(BR 1202-0725), 0. 20-0.23 Jy kms" 1 (PSS J2322+1944), 
and 0.13-0.20 Jy kms" 1 (APM 08279+5255) and are in 
good agreement with our observations for all sources. 
As described in Sect. 3.2.2, the observed 12 CO( J=l->0) 
linewidth in our new observations and in the VLA spec- 
trum (Carilli et al. 2002a) of PSS J2322+1944 is lower 
than that of the higher-J transitions, thus the integrated 
model fluxes are a bit higher than our result. If we com- 
pare the peak fluxes only, we obtain a good agreement. 

Our observations and models are in agreement with the 
assertion that all observed 12 CO( J=l— >0) flux density is 
associated with the highly excited molecular gas seen in 
the high- J CO lines. We thus find no evidence for an ad- 
ditional luminous, more extended low surface brightness 
gas component surrounding the central region of our tar- 
get QSOs, in contrast to what has been suggested previ- 
ously for APM 08279+5255 (Papadopoulos et al. 2001). 
Given the accuracy of our measurements, we conclude 
that at most 20-30% of the 12 CO(J=1^0) luminosity 
may be associated with such a diffuse component. We 
note however that, if the L' co to M ga8 (H2) conversion 
factor ('a') for a faint extended component were higher 
(e.g. Galactic), a higher H2 mass may be hidden in such 
an extended component. 

4.2. Correlations of high and low redshift galaxies 

CO results for our three high-redshift IR-luminous 
QSOs are summarized in Tab. ^ and [5] Our sample 
consists of all three high-z QSOs for which the CO 
ground-state transition has been detected to date, and 
covers ~20% of the CO-detected z > 2 quasars. As 
discussed above, the 12 CO(J=1^0) transition provides 
the best information of the total amount of molecu- 
lar gas in a system, quantified by -^co(i-o)' ^ or mos t 
QSOs/galaxies detected in CO at high z, 12 CO( J=3^2) 
and 12 CO(J=4^3) are the lowest transitions that have 
been detected to date. Based on our results in the pre- 
vious paragraph, however, we can now estimate their 
^co(i-o) k v assuming constant brightness temperature 
for all transitions from 12 CO(J=1^0) to 12 CO( J=4-^3) 
(i.e., L' co is the same for those transitions). This as- 
sumption is also in agreement with observations towards 
the radio galaxy 4C60.07 (Greve et al. 2004). However, it 
is important to note that all high-z sources were selected 
via higher-J CO transitions, which could, in principle, 
introduce a bias towards highly excitated starburst envi- 
ronments. 

Following Sanders et al. (1991), Gao & Solomon (2004) 
found a non-linear relation between the logarithms of the 
FIR luminosity Lfir and the 12 CO(J=1^0) line lumi- 
nosity L'qq for a sample of local spiral galaxies, luminous, 
and ultra-luminous infrared galaxies. Their sample con- 
sists mostly of galaxies whose FIR luminosity is powered 
by star formation only. Our 12 CO( J=l— >0) observations 
and LVG models suggest that L' co can be estimated for 
all CO-detected high- 2: sources with some degree of con- 
fidence, even if only observations of higher-J CO transi- 
tions exist for most of these sources. To put this result 
into context, we now aim to discuss the J^o^FiR re ~ 
lation and its implications for some selected samples of 
galaxies at low and high redshift. 

Figure [5] shows the relationship between log(JFm) and 
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\og(L' co ) for our three targets, all other high-z CO de- 
tections (except TN J0924-2201, which does not have a 
measured FIR luminosity), the z < 0.2 Palomar-Green 
(PG) QSOs from Alloin et al. (1992), Evans et al. (2001), 
and Scoville et al. (2003) including PDS456 (Yun et 
al. 2004), local (z < 0.3) ULIRGs from Solomon et al. 
(1997), and the local galaxy sample (z < 0.1) of Gao 
& Solomon (2004). All FIR luminosities for the high- 
z sources are re-derived as described by Carilli et al. 
(2005) unless stated otherwise. Lensing magnification 
factors were taken into account. In addition, data for 
the Milky Way disk (Fixsen et al. 1999) and an (extrap- 
olated) relation for Galactic molecular clouds (Mooney 
& Solomon 1988) are given for comparison. The solid 
line is a straight-line least squares fit to the Gao & 
Solomon (2004) sample, corresponding to log(LpiR) = 
(1.26 ± 0.08) x log(-Lco) — 0.81 (the power law index is 
1.25 ± 0.08 in the original publication). Fitting the high- 
z sources, the Gao & Solomon (2004) data, the Solomon 
et al. (1997) ULIRGs, and the PG QSOs together, we find 
the relation log(L Fm ) = (1.39 ±0.05) x log(Z C o) - 1.76. 
Both relations are clearly non-linear, as their power-law 
indexes are significantly larger than unity (see also the 
corresponding discussion in Gao & Solomon 2004). We 
recover a Schmidt-Kennicutt law (power law index of 
1.4, Kennicutt 1998a, b) for the larger, more heteroge- 
nous sample. More importantly, the high-z sources seem 
to follow the same slope (within the errors) as seen at 
low z, even though most high-z sources harbor central 
AGNs. In addition, the high-z QSOs and high-z sub- 
millimeter and radio galaxies (likely lacking a luminous 
AGN) statistically occupy the same area in the plot. The 
same behaviour is found for the PG QSOs in compari- 
son to the low-z spirals/LIRGs/ULIRGs. It is in general 
not clear that the FIR emission comes mostly from star- 
formation, i.e. that it is not strongly biased by active 
nuclei or the interstellar radiation field. E.g., SED mod- 
eling of the z = 2.6 quasar H1413+117 (the Cloverleaf) 
reveals that only ^20% of its FIR luminosity is pow- 
ered by star formation (Solomon et al. 2003), but it is 
not an outlier in the L'qq-Lyir. diagram. In this con- 
text, the elevated L Fm /L' co ratio in APM 08279+5255 
may be explained by a combination of differential lens- 
ing and a relatively high contribution of the AGN to 
^fir (see Rowan-Robinson 2000 for models of the IR 
SED). Outliers like APM 08279+5255 are also found at 
low z (e.g. the QSO PDS456 and the warm ULIRG 
IRAS 08572+3915, see Fig. EJ). As the relative number 
of such outliers is very small, this may also be explained 
by a brief FIR-bright AGN phase (Yun et al. 2004). 
However, it is likely that the dominant energy source 
in most ULIRGs is an extreme starburst rather than 
heating by a dust-enshrouded AGN (e.g., Solomon et 
al. 1997; Downes & Solomon 1998). This picture is sup- 
ported by the finding that ULIRGs harbor large quan- 
tities of dense molecular gas, which is more intimately 
involved with star formation than the major fraction of 
the mostly diffuse CO (Gao & Solomon 2004). It has to 
be kept in mind that the high-z sources are highly se- 
lected, and probably fulfill the Malmquist bias (i.e., the 
apparent increase in the average L' C q and Lfir towards 
high redshift is probably a consequence of the flux limita- 
tion in the sample, e.g. Sandage 1994). As an additional 
consequence of flux limitation in the sample, lensing in- 



fluences the results (i.e., stronger lensing magnification 
allows to probe deeper, and therefore the intrinsic CO 
luminosity drops). 

We can also set an approximate lower limit on the du- 
ration of the intense starburst phase. Considering the 
total H2 masses and star-formation rates in Table [21 the 
depletion timescale of the molecular gas is of order 10 7 yr 
for all three QSOs under the assumption of a constant 
SFR and 100% star-forming efficiency. This implies that 
the starburst itself can be relatively short-lived, (and 
compact, as the dynamical time must be less than the 
starburst lifetime) unless the molecular gas in which the 
star-formation occurs can be re-supplied on timescales 
of-10 7 yr. 

5. SUMMARY 

We have detected 12 CO(J=1^0) emission in three 
QSOs at redshifts 3.9 < z < 4.7 with the NRAO 
GBT and the MPIfR Effelsberg telescope. From our 
analysis of the resulting spectra of BR 1202-0725, 
PSS 2322+1944, and APM 08279+5255, we obtain the 
following results: 

1. We derived lensing-corrected 12 CO( J=l— >0) line 
luminosities of 1.0 x 10 n K kms^pc 2 for BR 1202- 
0725, 4.2 x 10 10 K kms-'pc 2 for PSS 2322+1944, and 
1.4 x 10 10 K kms-ipc 2 for APM 08279+5255. These re- 
sults are in good agreement (within a factor of 2) with 
previous estimates of the total CO luminosities based on 
the higher-J CO transitions, consistently providing very 
large M gas (H 2 ) of > 10 10 M Q . 

2. 12 CO(J=1^0) fluxes predicted by one-component 
LVG models are in good agreement with our observa- 
tions. Considering our modeling results, the CO emis- 
sion appears to be close to thermalized up to the 4— >3 
transition in all cases. Thus, our observations show no 
indication of a luminous extended, low surface bright- 
ness molecular gas component in any of the high-redshift 
QSOs in our sample (cf. Papadopoulos et al. 2001). In 
fact, all CO transitions are described very well by a single 
gas component where all molecular gas is concentrated in 
a compact circumnuclear region. If such extended com- 
ponents were to exist, our observations and models limit 
their contribution to the overall luminosity to at most 
20-30%. 

3. log(XpiR) and \og(L' co ) appear to be correlated for 
low-redshift galaxies over orders of magnitude in Lfir in- 
cluding (U)LIRGs, and the significantly brighter sources 
found at high-z appear to follow the same general trend 
(see also Solomon & Vanden Bout 2005). In particu- 
lar, we find that the correlation shows no significant dif- 
ference between QSOs and systems without a luminous 
AGN. 

The observations presented herein demonstrate the fea- 
sibility of detecting high-z CO with 100 m single-dish 
radio telescopes, and highlight the physical implications 
of observing the ground-state transition of this molecule 
towards massive galaxies at redshifts greater than 4. 
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Fig. 1.— GBT (top) and Effelsberg (middle) spectra of 12 CO(J=1^0) emission from BR 1202-0725. The bottom panel shows the 
combined spectrum. The GBT spectrum has been smoothed to a resolution of 2.00MHz (30 kms" 1 ). The rms per channel is ^75 fijy. The 
Effelsberg spectrum has been smoothed to a resolution of 5.97 MHz (88 kms" 1 ). The rms per channel is ~155 £ijy. For illustration, the 
inset shows the full spectral range of the Effelsberg spectrum that was used for spectral baseline fitting (width: ~380MHz). The combined 
spectrum has been smoothed to a resolution of 2.00 MHz (30 kms" 1 ). The rms per channel is ~70 fijy. The thin black lines show Gaussian 
fits to the data (see Tab. HI . 
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Fig. 2.— GBT (top) and Effelsberg (middle) spectra of 12 CO(J=1^0) emission from PSS J2322+1944. The bottom panel shows the 
combined spectrum. The GBT spectrum has been smoothed to a resolution of 1.80 MHz (24 kms -1 ). The rms per channel is ~140/^Jy. 
The Effelsberg spectrum has been smoothed to a resolution of 3.60 MHz (48 kms -1 ). The rms per channel is ~380 fijy. For illustration, 
the inset shows the full spectral range of the Effelsberg spectrum that was used for spectral baseline fitting (width: ~180MHz). The 
combined spectrum has been smoothed to a resolution of 1.80 MHz (24 kms -1 ). The rms per channel is ~125^tJy. The thin black lines 
show Gaussian fits to the data (see Tab. HI . 
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Fig. 3.— GBT spectrum of 12 CO(J=1^0) emission from APM 08279+5255. The spectrum has been smoothed to a resolution of 
5.86 MHz (75 kms -1 ). The rms per channel is ~65 pjy. The thin black line shows a Gaussian fit to the data (see Tab. HI. 
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Fig. 4. — CO ladder and LVG models (based on all J > 1 transitions) for BR 1202-0725. The inset shows a zoomed— in version of the 
J = 2, J = 1 region (as indicated by the dashed box). The filled triangle is the 12 CO( J=l— >0) data point based on the combined spectrum. 
Data for the higher— J CO transitions (open triangles) are taken from the literature (see Table [3J. For comparison, we also show data for 
the inner disk of the Milky Way (open squares, Fixsen et al. 1999), normalized to the 12 CO(J=l— >0) flux of BR 1202-0725. The kinetic 
temperature T^ n and gas density p gas (H2) are treated as free parameters in this study. Three representative models are shown: Model 1 
(solid line) assumes T^ n = 60 K and p ga s(H2) = 10 41 cm -3 , and gives the overall best fit to all transitions. Model 2 (dashed line) assumes 
T kin = 120K and p g as(H 2 ) = 10 3 7 cm- 3 , while model 3 (dotted line) assumes T kin = 30K and p ga s(H 2 ) = 10 4 6 cm- 3 . 
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Fig. 5. — Comparison of velocity— integrated CO line luminosity with FIR luminosity for a sample of low— z spiral and starburst galaxies 
from Gao k, Solomon (2004), the ULIRGs from Solomon et al. (1997), z < 0.2 Palomar-Green (PG) QSOs from Alloin et al. (1992), 
Evans et al. (2001), and Scoville et al. (2003), extrapolated Galactic molecular clouds (GMCs) from Mooncy & Solomon (1988), the Milky 
Way (Fixscn et al. 1999), and for high— z sub-millimeter galaxies (SMM), radio galaxies, and QSOs from the literature (see Solomon & 
Vandcn Bout 2005) with respect to our new results. All FIR luminosities are re-derived as described in Carilli et al. (2005). All data are 
corrected for gravitational lensing. The solid line is a straight— line least squares fit to the Gao & Solomon (2004) sample, corresponding to 
log(L F iR.) = (1.26 ±0.08) X log(Lco) -0.81. The dotted line shows the fit to the GMCs from Mooney & Solomon (1988), corresponding to 
a power law index of 0.96 ± 0.08. The dashed line is a fit to all high— z sources, the Gao & Solomon (2004) data, the Solomon et al. (1997) 
ULIRGs, and the PG QSOs, corresponding to log(L F m) = (1-39 ± 0.05) x log(L C o) - 1-76. 



TABLE 1 

Observed 12 CO (J— 1— >0) line parameters 



Source 


2 


^obs 

[GHz] 


[mjy] 


AVpwHM 

[kms- 1 ] 


[Jy kms -1 ] 


Telescope 


BR 1202-0725 


4.6932 ± 0.0004 
4.6956 ± 0.0012 
4.6949 ± 0.0003 


20.2450 
20.2388 
20.2411 


0.36 ± 0.03 
0.40 ± 0.11 
0.34 ± 0.03 


329 ± 36 
522 ± 146 
333 ± 30 


0.124 ± 0.012 

0.22 ± 0.06 
0.120 ± 0.010 


GBT 
Effclsberg 
(combined) 


PSSJ2322+1944 


4.1173 ± 0.0003 
4.1184 ± 0.0008 
4.1179 ± 0.0002 


22.5258 
22.5211 
22.5231 


0.77 ± 0.07 
1.29 ± 0.26 
0.77 ± 0.07 


202 ± 17 
184 ± 46 
190 ± 14 


0.165 ± 0.014 

0.25 ± 0.06 
0.155 ± 0.013 


GBT 

Effclsberg 
(combined) 


APM 08279+5255 


3.9122 ± 0.0007 


23.4663 


0.26 ± 0.04 


556 ± 55 


0.152 ± 0.020 


GBT 



Note. — All given uncertainties are formal (correlated) errors from the Gaussian fit. The z error was derived from a full 3 parameter fit. In the 
error estimates for S u and AVfwhm, the correlation with the error in z is neglected, as the relative uncertainty in z is by ~3 orders of magnitudes 
less as that of the other two quantities. The error estimate for / takes the correlation of the errors of S v and AVfwhm into account. 



a A T£(K)/S'(Jy) conversion factor of 1.5 was assumed for the GBT. 
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TABLE 2 

Luminosities, gas masses and Star-formation rates (SFRs). 



Source 


[Gpc] 




iFIR a 

[10 12 L ] 


J 1 a 
^CO(l-O) 
[10 1() K kms-ipc 2 ] 


^fir/^CO(i-0) 
[L q /K kms-ipc 2 ] 


M gas (H 2 ) b 
[10 10 M Q ] 


SFR C 
[H^MQyr- 1 ] 


BR 1202-0725 


44.2 


l d 


60 g /60 


10.1/10.1 


596 


8.1 


9.0 


PSS J2322+1944 


37.8 


2.5° 


23 h /9.2 


10.5/4.2 


219 


3.4 


1.4 


APM 08279+5255 


35.6 


7 f 


200728.6 


9.6/1.4 


2090 


1.1 


4.3 



a Apparcnt luminosities (not corrected for lensing)/intrinsic luminosities (lensing-correctcd) . 

b Assuming a conversion factor of a — 0.8Mq /K kms _1 pc 2 from ^co(l-o) ^° -^gas(H2) as appropriate for ULIRGs (see Downes & Solomon 
1998). 

c Assuming a Schmidt-Kcnnicutt law (Kennicutt 1998a, b): SFR[Mq yr~ 1 ]— 1.5 X 10~ 10 Lfir[L©], i.e., (5mf^sb — 1-5 following the notation of 
Omont et al. (2001): <5mf describes the dependence on the mass function of the stellar population, <5sb gives the fraction of Lrir that is actually 
powered by the starburst and not the AGN. 

d Carilli ct al. (2002b); 

°Carilli ct al. (2003); 

f Lewis et al. (2002); 

SCarilli ct al. (2005); 

h Cox et al. (2002); 

'Bcclcn ct al. (2006). 



TABLE 3 

CO DETECTIONS IN THE QUASAR HOSTS OF OUR TARGETS IN THE LITERATURE. 



Source 




2 


Transition 


[mjy] 


I 

[Jy kms -1 ] 


AVfwhm 
[kms- 1 ] 


ref. 


BR 1202-0725 


N 


4.692 


2-1 


0.44 ± 0.07 


0.26 ± 0.05 




1 




S 


4.695 


2-1 


0.77 ± 0.10 


0.23 ± 0.04 




1 




NS 




4-3 


~5.1 


1.50 ± 0.3 


280 ± 30 


2 




N 


4.6916 


5-4 


~3.5 


1.3 ± 0.2 


350 


2 




S 


4.6947 


5-4 


~5.5 


1.1 ± 0.3 


190 


2 




NS 


4.695 


5-4 


9.3 ± 2.1 


2.7 ± 0.41 


220 ± 74 


3 




NS 


4.6915 


7-6 


~10.6 


3.1 ± 0.86 


~275 


2 


PSS J2322+1944 




4.1192 


1-0 


0.89 ± 0.22 


0.19 ± 0.08 


200 ± 70 


4 






4.1192 


2-1 


2.70 ± 0.24 


0.92 ± 0.30 


280 ± 42 


4,5 






4.1199 


4-3 


10.5 


4.21 ± 0.40 


375 ± 41 


6 






4.1199 


5- 4 

6- 5 

7- 6 


12 


3.74 ± 0.56 


273 ± 50 


6 
7 
7 


APM 08279+5255 




3.9 
3.9 


1- 

2- 1 




0.150 ± 0.045 a ; 0.22 ± 0.05 
~0.81 b 


575 


8,9,10 
8 






3.9114 


4-3 
6-5 


7.4 ± 1.0 


3.7 ± 0.5 


480 ± 35 


11 
7 






3.9109 


9-8 


17.9 ± 1.4 


9.1 ± 0.8 
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References. — [1] Carilli et al. (2002b); [2] Omont et al. (1996); [3] Ohta et al. (1996); [4] Carilli et al. (2002a); [5] Carilli ct al. (2003); [6] 
Cox ct al. (2002); [7] Weiss ct al. (2006), in prep.; [8] Papadopoulos et al. (2001); [9] Lewis et al. (2002); [10] Riechers ct al. (2006), in prep.; [11] 
Downcs et al. (1999) 

NOTE. — For BR 1202-0725, N and S indicate the northern and southern component, NS indicates an integrated measurement over both 
components. 

a Dcrivcd for central ~ l" only. 

b Assuming a velocity-averaged brightness temperature ratio between 12 CO(J— 2— >1) and 12 CO(J— 1^0) of 1.35 ± 0.55 (Papadopoulos et al. 
2001). 



